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We study the relationship between the differential conductance and the local density of states in tight- 
binding tunnel junctions where the junction' geometry can be varied between the point-contact and the 
planar-contact limits. The conductances are found to differ significantly in these two limiting cases. We 
also examine how the matrix element influences the tunneling characteristics and produces contrast in a 
simple model of scanning tunneling microscope (STM). Some implications regarding the interpretation of 
STM spectroscopic data in the cuprates are discussed. The calculations are carried out within the real-space 
Keldysh formalism. 
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I. INTRODUCTION 

The spectrum of one-particle excitations is a rich source 
of information in the study of condensed-matter systems.^ 
Among the various tools available to measure the low-energy 
excitations, tunneling has become a prominent technique 
for electronic materials, due to its high energy resolution 
and the possibility to set up local probes.'^ In a tunneling 
experiment, electrons are extracted from or injected into 
a sample through a classically forbidden barrier, thereby 
probing the spectrum of excitations below and above the 
Fermi energy. The practical realization requires to keep 
the sample and the probing electrode sufficiently isolated 
that they do not significantly interact, and sufficiently close 
that a measurable number of electrons can tunnel from one 
to the other. This has been achieved in a variety of ways, 
including the planar junction with a thin insulating layer^ 
and the vacuum junctions obtained, e.g., by breaking the 
sample^'^ or using a scanning tunneling microscope (STM) 
tip.*""^ The precise relationship between the current -voltage 
characteristics measured at the tunnel junction and the 
spectrum of excitations has been a matter of continuous 
research since half a century. ^"^^ 

The problem may be divided into three distinct questions: 
(i) is the tunneling-Hamiltonian formalism^ ^ appropriate, 
and if yes, how should the tunneling Hamiltonian be de- 
fined? More elaborate approaches have suggested that the 
tunneling Hamiltonian can be a good approximation,^^'^''^° 
but the proper definition of the tunneling matrix element 
remains a difficult problem lacking a general solution, that 
would remain valid for interacting systems and out-of- 
equilibrium conditions, (ii) Provided that the tunneling- 
Hamiltonian gives a good description of the system and 
that the tunneling matrix elements are known, how does 
the current relate to the excitation spectrum in the elec- 
trodes? In particular, is it possible to measure the density 
of states (DOS) or the local DOS (LDOS) when using local 
probes? This question is crucial for the interpretation of ex- 
perimental data, but has not been thoroughly investigated, 
(iii) Although surfaces are of great interest by themselves, 
tunneling spectroscopy is often conducted with the aim of 
addressing bulk properties. Therefore, when it comes to the 
interpretation of tunneling spectra, the question invariably 



arises, whether the measurements performed by connecting 
two surfaces are representative of the bulk materials. We 
will not address (i) and (iii), but turn our attention to (ii), 
which is a well-posed theoretical problem. In order to cope 
with (i), we start right away from a tight-binding Hamilto- 
nian containing two electrodes and a tunneling term. In 
order to avoid (iii), we shall study systems in which the 
distinction between the bulk and the surface is irrelevant, 
at least for the electrode representing the sample. 



The goal of the present study is to explore, using mini- 
mal models, how the junction' geometry and the specifics 
of the tunneling term influence the current-voltage charac- 
teristics I(y). We are particularly interested in comparing 
the differential conductance cr(V) = dl/dV and the DOS, 
since these two quantities are commonly believed to be 
roughly proportional to each other. The DOS and, more gen- 
erally, the one-electron excitations, are very conveniently 
described by means of the Green's function, which is an 
energy-dependent nonlocal quantity G{x,x', co) in the real- 
space representation.^ The LDOS, in particular, is directly 
related to the diagonal elements of G. The real-space formu- 
lation is best suited for our purposes, as we shall investigate 
geometrical effects and consider systems that break transla- 
tion invariance. The exact current I(y) can be expressed in 
terms of the real-space Green's functions using the Keldysh 
nonequilibrium formalism.^^ In Sec. II, we recall this for- 
malism, which leads to a compact formula for the current 
[Eqs. (2) and (3)] involving only the Green's functions of 
the electrodes and the tunneling amplitudes. ^■^'^^ In Sees. Ill 
and IV, we use this formula to compute the /(V) characteris- 
tics of two simple tight-binding models designed in order to 
highlight the role of geometry and dimensionality on the re- 
lationship between cr(V) and the DOS. These models show, 
in particular, that Harrison's cancellation argument,^ which 
explains the ohmic behavior of planar junctions between 
simple metals, has a relatively narrow range of validity. In 
Sec. V, we introduce a toy model of STM and discuss the role 
of the matrix element in the imaging process and the local 
spectroscopy. Finally, our main results and their implications 
are summarized and discussed in Sec. VI. 
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II. MODEL AND FORMALISM 

The hallmark of vacuum tunnel junctions, when com- 
pared to other electrical contacts, is that the two electrodes 
remain electrically well isolated when the junction is formed. 
This is due to the presence of a high barrier strongly reduc- 
ing the wave-function overlap. In this situation, a good 
starting point is to represent the junction by a tunnel Hamil- 
tonian that describes single-electron transfers across the 
barrier. The electrodes, on the other hand, are assumed to 
be in thermal equilibrium, and the Hamiltonian describing 
them to remain unchanged in the presence of the junction, 
apart from the electrical bias applied to one of the contacts. 
This assumption can only be valid if the potential drop of 
the biased junction takes place entirely in the region of the 
barrier The extent to which these various assumptions apply 
to realistic experimental setups is obviously a difficult one, 
and we do not aim to address it here. We consider instead a 
class of ideal tight-binding junctions where these assump- 
tions are build in from the start. The junction's differential 
conductance can then be evaluated exactly without further 
approximation, and the relation between the conductance 
and the LDOS can be inspected. In this section, we first 
review the general formula for the conductance in terms of 
the Green's functions of the two electrodes and the tunnel- 
ing amplitudes, and then consider a few remarkable special 
situations where the expression of the conductance can be 
simplified further Although we restrict our applications to 
the case of noninteracting single-band metals, many of the 
quoted expressions apply equally to multiband interacting 
systems, as long as the correlations induced by the tunnel- 
ing term are negligible in comparison to the correlations 
intrinsic to the electrodes. 



A. General formula for the differential conductance 

Consider two electrodes, left (L) and right (R), described 
by the lattice tight-binding Hamiltonians and respec- 
tively. The electrodes being initially decoupled means that 
these Hamiltonians commute: [J^^,^] = 0. We denote 
the lattice sites in L and R by I and r, and the correspond- 
ing electron field operators by i/'^CO and ip'Xr). The spin 
coordinate of the electrons is omitted throughout. The two 
electrodes are coupled by a tunneling term 



ni,r)^p\l)i^ir) + H.c., 



(1) 



so that the system's Hamiltonian is = M'l + M'^ + M'j . 
The tunneling amplitudes T[l,r) are symmetric, and we 
also assume that they are real. When a bias is applied to 
L, a steady-state current / is established. The calculation of 
the current is outlined in Appendix A and gives 



dwReTrTG^Cw). 



(2) 



T is the matrix with elements T{1, r), and is the Keldysh 
matrix Green's function defined in the time domain by 



G^{r,l,t) = -i{[i/j(r,t),i/j'^(;,0)]) with {•••) the thermo- 
dynamic average with respect to M', the fully formed junc- 
tion. The trace is over the sites I of L. Hence the current 
sums all closed paths from a point I to itself with excursion 
in R through the barrier For convenience, it is desirable 
to express G^(a)) in terms of the Green's functions of the 
isolated electrodes. At this step, we must implement the as- 
sumptions that the bias-induced potential drop takes place 
in the barrier region, and that the two electrodes are in 
thermal equilibrium at the chemical potentials /i^ and /i^. 
This gives (see Appendix A) 

TG"^ = (a - TG+TG+T^ [(1 - 2/i)rG+r(G+ - G-) 

+ (1 - 2/R)r(G+ - G-)rG-] (1 - TG-TGir'- (3) 

G^^(a)) are the retarded (-I-) and advanced (— ) matrix 
Green's functions in L and R. In the time domain, they read 

G±ai,i2,t) = Ti0(±t)(W(ii,o,i/''a2,o)]+) with [•,•]+ 

being the anticommutator, and, similarly, for G^. In the 
frequency domain, we have G^^(£o) = [G^jjCco)]'". fi^R(.co) 
is the Fermi distribution measured from /^^ jj. Equations (2) 
and (3) allow to evaluate the current-voltage characteristics 
/(V) as well as the tunneling conductance dl/dV once the 
electrodes and the tunneling amplitudes are specified. The 
size of the matrices to be considered is given by the number 
of sites connected by nonvanishing tunneling matrix ele- 
ments. The propagators Gf^[co) must be calculated for the 
semiinfinite electrodes, but they are only needed among this 
same subset of sites for calculating the current. It should be 
noted that while Eq. (2) is general, Eq. (3) is not exact if 
there are interactions in the electrodes. For noninteracting 
electrodes, Eqs. (2) and (3) reduce to Eq. (37) of Ref. 22 
(see Appendix B). The derivation of Ref. 22 explicitly relies 
on the assumption of independent electrons, and seems diffi- 
cult to generalize. An advantage of the Keldysh formulation 
is that interactions can be easily included, at least formally. 
In the case of interacting electrodes, it can be, for example, 
shown that the current obtained by inserting the interact- 
ing Green's functions in Eq. (3) correctly accounts for the 
correlations intrinsic to the electrodes, but disregards the 
correlations induced by the tunneling term. Within their 
domain of validity, Eqs. (2) and (3) are exact at all orders 
mT[l,r) and V and for arbitrary temperatures. 

In vacuum tunneling experiments, one generally uses a 
probe to measure a sample. From now on, we shall con- 
sider L as the probe and R as the sample. It is convenient 
to measure energies from the chemical potential of the 
sample: we thus take = and fi^ = eV. As a result, 
fn^co) = [expC^w) -I- = /((u) with j3 the inverse tem- 
perature, and fi^co) = /(w — eV). With this convention, 
G^(£o) to be used in Eq. (3) is the same as that for the 
isolated sample R, while G^(£o) should be the Green's func- 
tion in L with the origin of energies shifted by eV, i.e., 
gJ(£o) = G^g(£D — eV) with G"^^ the Green's function in 
the absence of bias. Therefore two kinds of terms appear 
when taking the bias derivative dl/dV: those from the bias 
dependence of fi, and those from the bias dependence of 
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FIG. 1. Generic tight-binding tunnel junction (left), and particular examples with the corresponding tunneling matrix element. The black 
dots represent discrete points of an arbitrary tight-binding lattice. 



Gf. For L to be a useful probe, its properties must vary 
slowly with energy, so that the latter terms are expected to 
be small compared to the former We consequently split the 
differential conductance and define 



dl 
dV 



a(y) + 5a(y). 



Retaining only the bias dependence of in Eq. (3), the 
dominant contribution to the conductance is 



dec [-/'(w - eV)]2ReTrrG+r(G7 - G+) 
xiTL-TG-TG-T^iTL-TG+TG+T^. (5) 



The o) dependence of G^^j is implicit, /' is the derivative of 
/, and we have used the cyclic property of the trace. Equa- 
tion (5) is a convenient starting point for performing analyti- 
cal calculations: at zero temperature, — /'( w — eV) becomes 
5(w — eV), and the energy integral drops. o"(V) then only 
depends on G^(eV) = Gf^^^O) and G^(eV). It is also seen 
that the main effect of a finite temperature on o-(V), apart 
from small contributions related to the bias dependence of 
G^, is to thermally broaden the zero-temperature curve. We 
have checked on simple models that <5-(V) gives the same 
conductance as the Landauer-Biittiker formula^^ at zero 
temperature, provided that the bias eV in L is taken into ac- 
count when calculating the scattering matrix. Hence 5cr(V) 
corresponds to the difference between the differential and 
the ballistic conductances. 



B. Particular cases 

In this section, we present more explicit expressions for 
the conductance, in a few limiting situations depicted in 
Fig. 1. While for a generic junction the matrix structure 
in Eq. (5) cannot be simplified further, in these cases the 
conductance can be reduced to a scalar form. The first 
case is the "point contact", by which we mean a junc- 
tion such that the tunneling is only possible from a sin- 
gle site Iq in L to a single site Tq in R. The tunneling 
amplitude is r(i,r) = t5;,^5rr„. In reciprocal space, the 
matrix element is therefore independent of momentum, 
iTj-fc/p = t^. The latter assumption is quite common — 
although often implicit — in theoretical studies of tunnel 



junctions. The fact that it is a drastic simplification ap- 
pears somewhat more clearly when it is formulated in real 
space. The locality of the tunneling amplitude implies that 
all matrices in Eq. (5) become scalars. Using the defini- 
tion of the LDOS, e.g., iVi(Io,w) = --ImG+(Io, Iq, w) = 



(4) ^ [G^ (io, Iq, (o) - G+(Io, lo, w)], one readily obtains 
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C7(V)=- dco[-fXco-eV)]NL[lo,co) 
h 

^ |l-t2G+ao,io,w)G+(ro,ro,w)|2 

(point contact). (6) 

At zero temperature and in the weak tunneling regime 
t ^ 0, we thus have an exact proportionality between 
o-(V) and the sample LDOS at the point Vq-. o"(y) oc 
iVi_o(io, 0)JVR(ro, eV) + ^(t^). N.^oik, 0) = N^Ho, eV) is the 
LDOS at the Fermi energy and at the point ig in L. This 
proportionality is lost as t increases and becomes compa- 
rable to the typical energy scales in L and R, so that the 
product t^G^G'^ in the denominator becomes of order unity. 
In the class of junctions that we are considering, the weak 
point-contact limit is the only case where there can be a 
strict proportionality between the tunneling conductance 
and the LDOS.^^ 

A second case of interest is when the tunneling is possible 
from a single point Iq in L to several points r in R. This 
can be considered as a crude model for an STM with some 
dispersion of the electrons tunneling from the tip. The 
matrix element reads r([, r) = 5;,^t(r). It is still possible to 
invert the matrices analytically in this case, and one obtains 



dcD[-f'{co-eV)]N^{lQ,CD) 



(271)2^ t(ri)(-^)lm[G+(ri, r^, oj)]t{r^) 



1 - G+([o, lo, oj) 2 t(ri)G+(ri, r^, w)t(r2) 



("STM" junction). (7) 

When compared to Eq. (6) — that is recovered by taking 
t(r) = t5rr„ — Eq. (7) includes all excursions within R 
among sites that are connected back to Iq. This leads to 
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interference between the various paths and to the loss of a 
simple proportionality between u(y) and the LDOS in R. 

The planar analog of the point contact is a junction 
such that tunneling is only possible between sites facing 
each other on the two surfaces. The matrix element is 
r(i,r) = t5ij^5^^^^5^,.^^, where Zj^ and denote the spa- 
tial coordinates in the direction perpendicular to the junc- 
tion's plane, Zg and are the coordinates of the surfaces of 
L and R, and and ry denote the coordinates in the plane 
of the junction. In the momentum representation fcy, all 
matrices in Eq. (5) again become scalars, leading to 



C7(V) = - 

h 



h 

(27it)%(fc||,ro,ro,tL)) 

^ |1 - t2G+(fc||, Zo, Zo, w)G+(fc||, ro, r^, 

(local planar junction). (8) 

We have introduced the local spectral functions 
Ai(fc||,Zi,Z2,w) = ^[G-(fc||,Zi,Z2,w) - G+(fc||,Zi,Z2,w)], 
etc., whose diagonal part A^(fc||,Zo,Zo, w) has the physical 
meaning of a fcy -resolved LDOS on the surface of L. 
The crucial difference between Eqs. (6) and (8) is the 
conservation of in-plane momentum in the latter, which 
imposes to match states with the same fcy across the barrier. 
This restriction implies that the conductance per channel 
must be smaller in the planar than in the point-contact 
case, as we shall see below. Finally, we may consider a 
straightforward generalization of the local planar junction 
by including additional tunneling paths as depicted in the 
rightmost panel of Fig. 1. The matrix element becomes 
r(i,r) = 5ij^5j.^^^t{l\\ — ry), but since the translation 
invariance is preserved the structure of the matrices remains 
the same as in the previous situation: 



e 



dco [-fXcj - eV)] ^AiCfcy, Zo, Zo, CO-) 
[27rt(fcy)]^Aj;(fcy,ro,ro,to) 

|1 - t2(fcy)G+(fcy , Zo, Zo, w)G+(fcy, To, To, 

(Planar junction). 



(9) 



with t(fcy) the Fourier transform of t([y — ry). 

In the subsequent sections, we shall use Eqs. (2)-(4) to 
compute numerically the exact conductance, including the 
small correction 5 a due to the electronic structure in L. For 
the discussion and interpretation of the results, however, 
extensive use will be made of the simpler expressions (5)- 
(9). 



III. TWO-DIMENSIONAL T-SHAPED JUNCTION 

We have seen that in the point-contact limit and in the 
weak-tunneling regime the conductance is proportional to 
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FIG. 2. Tunnel junction between a semiinfinite rod and an infinite 
one-dimensional chain. For = 1, the junction is a T-shaped point 
contact while for n^^ = oo, it is a local planar junction according to 
the terminology of Fig. 1. 



the sample LDOS. This proportionality is lost when the cou- 
pling increases. It is also expected to be lost in the planar 
geometry due to the conservation of in-plane momentum. 
In order to investigate these various limits and the transi- 
tion between them, we consider the system represented in 
Fig. 2. The probe L is made of interconnected semiinfi- 
nite chains forming a flat rod. The Hamiltonian Jff^ involves 
only nearest-neighbor hopping with energy . The sample 
R is an infinite linear chain with nearest-neighbor hopping 
t^. We consider L and R to be half- filled for simplicity. Fi- 
nally, the tunneling term J^j connects neighboring sites at 
the interface with the hopping t. For n^, = 1, we have a point 
contact, while for = oo, we have a local planar junction 
according to the terminology of Fig. 1 . This is perhaps the 
simplest model to address the relation between the conduc- 
tance and the LDOS. The sample has a position-independent 
LDOS given by 



Nr{co) ■- 



-Re- 



^l-[W(2t«)]^ 



} . 



(10) 



This LDOS has square-root band edges divergences at ±2| tfil, 
that are easily identified if they appear in the differential 
conductance. We need the retarded and advanced Green's 
functions in both electrodes. In R, they can be given in 
closed form: 



GR(ri,r2,w): 



dk 



,ik(,ri-r2) 



2n CO- 2txC0sfc± iO+ 



1 ( co±iO+ 

2t- 



•-R 



2t 



R 



(11) 



0(Rez) (^z - _ 1 j 
-0(-Rez) ^z+^/z^-l 



\n\ 



|n| 



5 



3 






1 ' 1 


1 

1 = 1 - 

1 = 2 


\/f. 1 


1 , 1 





FIG. 3. Local density of states on the surface of the left system 
in Fig. 2. As increases, the LDOS in the central portion of the 
surface approaches the limiting value for = oo, while the LDOS 
close to I = 1 has a different shape. 



On the surface of L (x = 0), we may express the Green's 
function in terms of the single-electron wave functions 
^kqix,y) as 



2n CO — eV — Sl^ ib lO"^ 



(12) 



Indices Z^, = 1, . . . , run on the chains in the y direction, 
e^q are the single-electron energies with k (q) the momen- 
tum along X (y), and we have introduced the bias eV as 
discussed in the previous section. For open boundary condi- 
tions, ifkqix, 0) = (fkqix, ni + l) = (/JfcqCl, Z) = 0, the wave 
functions that solve are 



(^,,(x,Z)= [e'^^C-i) _ e"^a-.)] 



n,+l 



sin(qZ). (13) 



The energies are Sj^^ = 2ti(cosk + cosq) and the discrete 

q momenta are q = A = l,...,n,. With this, the k 

integral can be evaluated and the Green's function becomes 

, 2/ti ^ 

G (Zj, Z2, co) = — - > sm(qZi) smCqZj) 

"i + 1 g 

fco-eV±iO+ A 
xL cosq (14) 

L(z) = z- y/z + ly/z-1. 

The matrices Gfj, of Eqs. (11) and (14), together with the 
matrix T, which is just t x H in this case, must be substituted 
into Eqs. (2)-(4) for evaluating the differential conductance. 

The LDOS at the surface of L can be obtained from 
as iV^(Z, co) = (— l/7i)ImG^(Z,Z, co); it is neither energy nor 
position independent due to the open boundaries. Since 
we use L as a probe, some understanding of NiQ,co) is 
necessary before discussing the conductance. In Fig. 3, we 
display N^Q, co) for various values of n^. At = 1, there is 
only one point in contact with the sample R. The LDOS at 
that point is smooth and takes the form of a half-circle, 
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V^l- [W(2tJ]', (15) 




FIG. 4. Zero-bias differential conductance of the T-shaped junction 
depicted in Fig. 2 for = lOt^; and zero temperature. The inset 
shows the same data rescaled as indicated. 



with a half-bandwidth of 2| | . This contrasts with the LDOS 
far in the chain (x <C 0), which develops band-edge square- 
root divergences and approaches the form given by Eq. (10) 
with ti substituted for t^. The absence of divergences in 
JV^(Z, co) qualifies L as a convenient probe. This remains true 
for > 1, although additional structures develop in the 
LDOS. For = 2, the two sites on the surface have the same 
LDOS by inversion symmetry, with cusps at co = ±|t^|, and 
the half-bandwidth increases to 3|t£^|. For = 3, the LDOS 
on the central site (Z = 2) and on the borders (Z = 1, Z = 3) 
have the same bandwidth but different energy dependencies, 
and the same applies to all ^ 3. The general trend is that 
the LDOS becomes smoother and smoother as increases, 
and approaches the limiting value for = 00 with a half- 
bandwidth of 4|ti|, except at the borders Z = 1 and Z = 
where a different limiting behavior is found. With increasing 
Til, the characteristic energy of the structures in the LDOS 
decreases, but at the same time their amplitude is reduced. 
For L to be a good probe at any n^, the structures in N^Q, co) 
should not prevent the identification in the conductance of 
structures due to the sample R. This can be achieved by 
taking » t^. 

We have calculated numerically the conductance of the 
T-shaped junction of Fig. 2 as a function of n^,, t, and V for 
ti = lOtjj and zero temperature. The zero-bias conductance 
is presented in Fig. 4. In the tunneling regime where t <C 
ti, tjj, cr(0) increases as and is roughly proportional to 
the contact size as shown in the inset. Specifically, we 
have 



a(0): 



e2 2t' 



+ 0[t^) 



(16) 



L>-R\ 



with a[ni) a slow function of n^: a(l) = 1, a(2) = V3/2, 
and a^n^ » 1) ^ 1. a(2) is equal to the ratio of the zero- 
energy LDOS in L for n^, = 2 and = 1. Hence the reduced 
conductance for = 2 reflects the reduced zero-energy 
LDOS seen in Fig. 3. In Fig. 4, one sees that for = 1, 
the zero-bias conductance is bounded by one conductance 
quantum. This limitation is imposed by the probe, which 
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FIG. 5. (a) Differential conductance of the T-shaped junction 
depicted in Fig. 2 for = lOtj; = lOOt and zero temperature 
(black lines), compared with the sample LDOS [thick gray line, 
Eq. (10)]. The dashed gray curve for n,^ = oo represents Eq. (20) 
and does not contain the small correction 5a{V). (b) Function 
\M^i^\^ defined in Eq. (19) for the lowest discrete momentum q = 
n/ln^ + !)■ The function is unity for = 1, and approaches 
2n5{k — q) as —>■ oo, imposing the conservation of momentum 
at the planar junction. 



offers only one conducting channel. For > 1, the limit 
comes from the two outgoing channels of the sample, and 
the conductance has an upper limit of 2e^/h irrespective 
of n^. Finally, at large t, o-(O) drops as 8(e^/h)|titjj|/t^ 
irrespective of the contact size . 

In Fig. 5(a), we compare the bias-dependence of the dif- 
ferential conductance cr(V) with the sample LDOS A/^CeV) 
in the tunneling regime t <C ti, tj?- For n^, = 1 and = 2, 
the two quantities are almost undistinguishable, but for 
Til > 2, they start to deviate. Most noticeably, the signa- 
ture of the band-edge singularity in JVg(eV) is suppressed 
from the conductance, and the contact becomes ohmic as 



we approach the planar junction limit = oo. 



choice of parameters, the correction 5cj to the conductance 
remains small compared to a in the range of biases con- 
sidered. Therefore, we shall neglect 5a for the analysis of 
the results in Fig. 5. At = 1, Eq. (6) applies and we can 
deduce at zero temperature: 



(V): 



47Tt2 



[l + 7zty\t,\NR^eV)] 



2 ■ 



(17) 



This very simple expression results because the real part of 
the in-chain local propagator G^(0, 0, to) vanishes within 
the band, so that only the imaginary part, namely the DOS, 
remains in the denominator. Equation (17) confirms that 
in the point-contact limit, the differential conductance is 
proportional to the sample LDOS for t ^ 0. According 
to Eq. (17), the conductance should vanish for V > 2|tjj|, 
while the actual conductance turns slightly negative above 
2|tsl (not shown in the figures). The negative 5a(y) has a 
simple interpretation: when the bias reaches the band edge 
at 2|tx|, the current starts to decrease due to the downward 
curvature of the probe LDOS (see Fig. 3), resulting in nega- 
tive differential conductance. This effect cannot be captured 
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FIG. 6. Differential conductance of the T-shaped junction depicted 
in Fig. 2 with = 1 (left panel) and ji^^ = 8 (right panel) for 
= lOtj; and various t. The thick gray lines show the sample 
LDOS Eq. (10). 



by o-(V), which neglects the bias dependence of the Green's 
functions in the probe. 

For arbitrary n^, the conductance can be evaluated explic- 
itly (see Appendix C) as 



where 



NR[eV)Y,\sinq\\Mgt\^ + 0[t^), (18) 



sin^q {l — cos [qCn^ -I- 1)] cos [fcCn^, + 1)] } 



qk\ 



(n^ -I- l)(cosq — cosfc)^ 

(19) 

and cosfc = eV/(2tR). The function IM^^P deals with the 
matching of momenta at the junction. It is displayed in 
Fig. 5(b) as a function of k for the lowest q value. At = 1, 



With our there is only one q value, q = 7i/2, and |M, 



qk\ 



1 SO that we 

recover Eq. (17). For = oo, we have |Mqj.p = 2Ti5{k — q) 
expressing the conservation of momentum at the planar 
junction. In this limit, the q sum in Eq. (18) can be converted 
to an integral and yields n^lsiufcl. Since fc is fixed by V, 
we furthermore have n2,|sinfc| = ni/[27T|ts|JVj;(eV)] for 
|eV/(2fi)| < 1 and n^l sinfc| = otherwise. Hence the LDOS 
iVs(eV) in Eq. (18) is exactly canceled in the planar junction 
and we get ohmic behavior: 



-0(2|t^|-|eV|) + ^(t4). (20) 



This result can also be derived directly from Eq. (8). It 
shows that the cancellation of the LDOS singularity in the 
planar junction requires the exact conservation of momen- 
tum at the interface. Although the ohmic conductance of 
Eq. (20) seems to be an instance of Harrison's cancellation 
argument,' it is not. In Harrison's treatment of planar junc- 
tions between simple metals, there is a cancellation between 
the group-velocity and the DOS in the direction of the current. 
The DOS factors in the direction perpendicular to the cur- 
rent are also canceled, but for a different reason, namely, the 
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exponential suppression of large perpendicular momenta 
in the tunneling matrix element. In our case, the matrix 
element is momentum independent, and furthermore the 
sample has no extension along the x direction, so that both 
ingredients of Harrison's argument are irrelevant. 

To finish this section, we point out that, in addition to 
the kinematic constraint due to conservation of in-plane mo- 
mentum, the increase of the tunneling amplitude also leads 
to a suppression of the LDOS singularity in the conductance, 
as illustrated in Fig. 6. Even in the point-contact limit, the 
band-edge singularity is suppressed as t increases. This is 
obvious from Eq. (17) since a^"''='^\V) ~ l/iVjj(eV) at large 
t. We stress, however, that the validity of our treatment 
becomes questionable as t increases, since the assumption 
that the potential drop is confined to the barrier region will 
loose validity at some point. For example, if t = t^, there 
is actually no tunnel barrier, and it is not clear where the 
potential drop should take place. 




FIG. 7. Tunnel junction to a two-dimensional layer. The probe 
(left) is semi-infinite in the z direction and has a finite cross section 
of size rij^^n^. The sample (right) is an infinite two-dimensional 
plane. For n^^ = 1 we have a point contact while for = oo we 
have a local planar junction following the terminology of Fig. 1. 



IV 3D JUNCTION TO A 2D LAYER 

In the previous section, we saw that the conservation of in- 
plane momentum and/ or the lowering of the tunnel barrier 
lead to a suppression of the sample band edge singularity in 
the tunneling conductance. Here, we investigate the case 
of a singularity that is not at the band edge, but at zero en- 
ergy. In the system of Fig. 2, one could induce a square-root 
divergence at zero energy by adding next-nearest-neighbor 
hopping to Jt^. The system of Fig. 2 has academic rather 
than practical interest, however We consider instead the 
junction represented in Fig. 7, which shares some geometri- 
cal similarities with realistic systems, in particular the tunnel 
junctions involving quasi-two-dimensional (2D) layered ma- 
terials like the high-T^ superconductors. The probe L is a 
bar of lateral size n^, x rii extending from z = to z = — oo, 
and the sample R is an infinite plane. Nearest-neighbor sites 
in L (R) are connected by a hopping energy t^, (t^), while 
neighboring sites at the junction are connected by a tunnel- 
ing amplitude t. As before, we take L and R half- filled for 
simplicity. 

In the sample, the Green's function reads 



(21) 



with = 2tj5(cosfc^ + cosfc^,) and rj 2 the sites of the 2D 
lattice. The corresponding LDOS (— l/7r)ImG^(r, r, w) is 
translation invariant and presents a logarithmic van Hove 
singularity at zero energy: 

JVsCto) = :r4r-lK(l - [w/(4t«)]') 0(4|tBl - (22) 



K is the complete elliptic integral of the first kind. On the 
z = surface of the probe, the Green's function for open 
boundary conditions is a straightforward generalization of 



Eq. (14): 



Gf(ii,l2,to)= . "^^^^^y ^ sin(q^Zi^)sin(q^Z2;,) 



(Hi + ly 



X sin(q Zj )sin(q Z2y)L 



co-eV-E^± iO^ 
2t^ 



(23) 



with fig = 2t£(cosqj. + cosq^) and 2 the positions of the 

surface sites. The LDOS on the surface presents no 
divergence, but cusps whose energy scale and amplitude 
decrease as increases, like those in Fig. 3. Taking = 
lOtff is enough to maintain the probe-related conductance 
correction 5a small with respect to a. 

In order to evaluate accurately the Green's function in 
Eq. (21), we perform anal)^ically the momentum integration 
along one direction, and numerically along the other, using 
fast Fourier transforms with JV = 2^^ discrete momenta. The 
infinitesimal 0"*" is replaced by 1/iV, which ensures a high 
energy resolution. 

The resulting conductance is displayed in Fig. 8 for = 
lOtjj = lOOt and various probe sizes n^. At = 1, the 
conductance follows accurately the LDOS consistently with 
Eq. (6), since the junction is a point contact. Something 
surprising happens at =2: the conductance is completely 
suppressed at the band edges, while the signature of the 
van Hove singularity remains at the band center. This is 
due to an interference between two tunneling paths, which 
is destructive at the band edges. According to the general 
formula (2), there are three types of paths contributing to 



the current at lowest order in t for n, 



2. In the first 



type, the electron starts from one site on the surface of L, 
tunnels to the sample R and immediately tunnels back to 
the starting point. This path involves only local propagators 
in L and R. In the second type, the electron tunnels into 
R, propagates to a neighboring site, tunnels back to L, and 
closes the loop by returning to the initial site. It turns out 
that this path does not contribute because the advanced and 
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eVl\tj^\ 



FIG. 8. Differential conductance of the junction depicted in Fig. 7 
for tj^ = lOtji = loot and zero temperature (black lines), compared 
with the sample LDOS [thick gray line, Eq. (22)]. The dashed gray 
curve for = oo represents Eq. (25) and does not contain the 
small correction 5u[y). 

retarded propagators in L are identical at zero energy and 
cancel [see Eq. (5)] for nearest-neighbor sites. Finally, in 
the third path, the electron propagates to a next-nearest- 
neighbor site in R after tunneling, tunnels back, and closes 
the loop. This latter path contributes to the current with the 
same amplitude as the first one, but with the opposite sign 
because the local and next-nearest neighbor propagators in 
L are opposite at zero energy. Collecting all terms we find 

^(n.=2)(y) = _ ^ i^im [g+(0, 0, e V) - G+(l, 1, e V)] 

+ 0{t^) (24) 

with G+(0, 0, eV) and G+(l, 1, eV) being the local and next- 
nearest-neighbor propagators in R, respectively. At the band 
edges, corresponding to (0, 0) and (ti, tz) states, these two 
quantities have the same imaginary part, explaining the sup- 
pression of the conductance. The behavior of the conduc- 
tance is qualitatively similar for all sizes 5^ 2: remarkably, 
the signature of the van Hove singularity remains sharp in 
the conductance, even in the planar junction limit. This 
contrasts with the suppression of the square-root singularity 
in Fig. 5. For = oo, the conductance can be evaluated 
analytically starting from Eq. (8); we find 

-F^(t4). (25) 

Equation (25) is plotted as the dashed gray line in Fig. 8: in 
the planar junction the band edges are cut, but the van Hove 
singularity remains. This behavior is reversed if instead of 
increasing we keep = 1 and increase t (see Fig. 9) . In 
this situation, the conductance remains large close to the 
band edges, while the van Hove singularity gets suppressed. 
As discussed at the end of the previous section, however, the 
validity of the model becomes questionable as t increases. 




-4 -2 2 4 



FIG. 9. Zero-temperature differential conductance of the junction 
depicted in Fig. 7 in the point-contact limit (n^ = 1) for t^ = lOtj; 
and various values of t. 

We can conclude that Harrison's result^ does not apply in 
the case of Fig. 7 either: extrapolating the results of Ref. 9, 
one would predict an ohmic behavior for the planar junction, 
while the actual characteristics is strongly nonlinear, and the 
effective bandwidth deduced from the conductance curve is 
two times smaller than the true sample bandwidth. Further 
discussion of these issues is provided in Sec. VI. 



V MOMENTUM-DEPENDENCE OF THE MATRIX ELEMENT 

In this section, we discuss the most prominent effects as- 
sociated with a momentum dependence of the tunneling ma- 
trix element. For this purpose, we consider the system repre- 
sented in Fig. 10. The sample is an infinite two-dimensional 
layer like in the preceding section. It's Green's function and 
local DOS are given by Eqs. (21) and (22), respectively. The 
probe is a semiinfinite chain corresponding to the = 1 
limit in Sees. Ill and IV. The tunneling is not constrained 



t(r)/ 




> 






h 











FIG. 10. "STM-like" junction to a two-dimensional layer. The 
probe is a semiinfinite nearest-neighbor chain with hopping t^^, 
and the sample is an infinite two-dimensional plane with hopping 
tj;. Tunneling can take place from the extremity of the probe to 
various sites in the sample. 
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to be local, but can take place from the extremity of L to 
various sites r in R with an amplitude t(r). Figure 10 is a 
particular case of the "STM-like" junction sketched in Fig. 1. 
We will focus on describing how the matrix element t(r) 
can change the simple proportionality, which is realized in 
the local limit, between the conductance and the LDOS. 

The main contribution to the conductance at zero tem- 
perature can be deduced from Eq. (7). Making use of the 
translational invariance in the sample we find 



0.3 



0.2 



0.1 



h \tr 



1 + 



j:„\m\^G+ik,eV) 



(26) 



(eV- 



■e^ + iO+)"i with : 



: 2tj{(cos + 



where G+(fc,eV) : 
cos ky),A[i{k,eV) = 5[eV — e^) is the spectral function, t{k) 
is the Fourier transform of t(r), and N is the number of 
k points. We shall restrict to the tunneling regime t(r) <C 
t^, tjj, in which case we have the simple result 



■f^) = IT ?T ^ V I t(fc)|2A«(fc, eV) + 0[t^). (27) 




Hence the conductance is an average over sample states in 
reciprocal space, each state being weighted by |t(fc)|^. A 
formula analogous to Eq. (27) has been occasionally used to 
analyze tunneling spectra of the cuprates."^ We stress here, 
though, that the validity of Eq. (27) does not extend beyond 
the simple model of Fig. 10. In particular, this expression 
cannot capture interference effects taking place within the 
probe, such as those responsible for the suppression of the 
band-edge conductance in Fig. 8. 



FIG. 11. (a) Zero-temperature differential conductance of the 
junction depicted in Fig. 10, with a tunneling amplitude given by 
tg/ V2 for one pair of neighboring sites, and zero otherwise (see 
the inset) . The momentum dependence of the matrix element is 
indicated above the graph; the sample LDOS is displayed in gray 
for comparison, (b) Same as (a) with a tunneling amplitude given 
by to/2 for four sites around a plaquette, and zero otherwise, (c) 
Map of the tip-sample distance obtained in constant-current mode 
using the matrix element (28) with k = 5. The distance is largest 
(white) above the lattice sites, (d) Corresponding conductance map 
evaluated at eV = 2tjj. The conductance is largest (white) between 



The simplest situation is the local limit t(r) = to5„„, or ^j^^ ^^^-^^ ^-^^^ ^^^^1 parameters are t, = 10t« = lOOt, 



t(fc) = tf). In this case, the conductance reflects accurately 
the DOS as we have seen in the previous section; this is also 
clear from Eq. (27) since [l/N)J^,^AR(k,eV) = Nj^^eV). 
The second case of interest is when the "tip" is located in- 
between two sites of the sample, above a nearest-neighbor 
bond. To begin with, it is most natural to assume that 
the tunneling amplitude is the same for the two nearest- 
neighbor sites below the tip, and zero otherwise. The re- 
sulting differential conductance is shown in Fig. 11(a), and 
deviates significantly from the DOS: it is suppressed at the 
lower band edge and enhanced at the upper edge. This 
is easily understood considering that the matrix element 
t(fc) ~ cos(fc^/2) or cos(fcj,/2) depending on the bond ori- 
entation. In either case, the matrix element vanishes at 
(ti, n) (lower band edge) and is largest at F = (0,0) (up- 
per band edge). The asymmetric conductance can also be 
interpreted as the interference between two types of tun- 
neling paths in real space. The first path involves the local 
propagators in L and R, while the second involves the lo- 
cal propagator in L and the nearest-neighbor propagator in 
R; we thus find ct(V) oc -Im[G+(0, 0, eV) + G+(l, 0, eV)] 
at lowest order in Iq. For a (ti, n) state G^(0, 0, eV) and 
G+(l, 0, eV) cancel, while for a F state they add up. 

Let's now assume that the tip is located above the cen- 
ter of a plaquette. The tunneling amplitudes are equal 



for the four sample sites forming the plaquette, and they 
vanish elsewhere. The resulting conductance displayed in 
Fig. 1 1 (b) is strongly asymmetric, this time the signature of 
the van Hove singularity is completely suppressed, owing 
to the matrix element being zero both at (ti, tz) and (n;,0). 
Considering the various tunneling paths, we get in this case 
ct(V) oc -Im[G+(0, 0, eV) + 2G+(1, 0, eV) + G+(l, 1, eV)] . 

The two extreme cases of Figs. 11(a) and 11(b) suggest 
that by scanning a tip over a discrete lattice, one can observe 
large variations of the local tunneling conductance if the 
"tunneling cone" is sufficiently narrow, i.e., if only the sites 
closest to the tip contribute significantly to the current. As a 
generalization of these two extreme cases, we now consider 
a tunneling amplitude that is a continuous function of the 
tip position I, decaying exponentially with the distance to 
the tip: 



t(r): 



-K(|r-!|-d„) 



(28) 



do is a reference distance set equal to the sample lattice 
parameter, k controls the opening of the tunneling cone. 
We may relate k to an effective work function </> of the 
junction through the relation 4>/Ej! = {Kjk^)^ ~ (jc/ti)^. 
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(f) thus varies between 0.4 and 2.5 times £p for k ranging 
between two and five in units of the sample lattice. The 
vertical tip-sample separation d (out-of-plane coordinate 
of I) is adjusted in order to keep the total current fixed, 
in the fashion of a constant-current STM experiment. The 
reference current is set to the value calculated at eV = 2t^ 
when the tip sits on top of a lattice site. For k- = 2, the 
tunneling cone is wide and the current averages over sev- 
eral paths. As a result, the conductance is almost uniform 
when scanning the tip over the lattice. Furthermore, it is 
strongly asymmetric, almost completely suppressed at neg- 
ative bias and peaked at positive bias. Indeed, as k- ^ 0, 
we approach the limit t(fc) = tgS^. q where tunneling is pro- 
hibited except for the F state at the upper band edge. The 
situation is more interesting for jc = 5, where the tunnel- 
ing cone is narrower and the matrix element interpolates 
smoothly between the situations shown in Figs. 1 1 (a) and 
11(b). In order to maintain the current when moving away 
from a lattice site, the tip has to approach the sample. The 
resulting height field displayed in Fig. 1 1 (c) has a corruga- 
tion Ad/d K 14%. Figure 11(d) shows a conductance map 
calculated simultaneously at eV = 2tg. Remarkably, the 
regions near the plaquette centers appear with the largest 
conductance, because the total current is fixed, and at these 
locations the conductance increases with increasing bias 
[Fig. 11(b)]. This can be understood as follows: assume 
that the conductance varies as C7(V) oc V" with a position- 
dependent exponent a, then the constant-current condition 
/(Vq) = Iq implies that c7(Vo) oc 1 -I- a. Therefore it appears 
that in constant-current mode, a conductance map gives 
an image of the bias dependence of the local conductance, 
rather than a horizontal cut of the LDOS. A map of the con- 
ductance at fixed tip height, in contrast, provides an image 
similar to Fig. 11(c). 



VI. DISCUSSION 

The two models studied in Sees. Ill and IV show that there 
can be significant qualitative differences between the con- 
ductances measured at planar and point-contact junctions 
to the same given sample. The conservation of in-plane 
momentum at the planar junction leads in one case to a 
suppression of the band-edge singularity (Sec. Ill), in the 
other case to the complete suppression of the conductance 
at the band edge without affecting the singularity at the 
band-center (Sec. IV) . One can conclude that the differen- 
tial conductance is not in general proportional to the sample 
DOS at planar junctions, although features of the DOS may 
be present in the conductance, like e.g., the van Hove sin- 
gularity in Fig. 8. There seems to be no simple rule to 
determine a priori whether a singularity of the DOS will or 
will not show up in the conductance of a planar junction. 
In contrast, at the point contact, a simple proportionality is 
found between the conductance and the sample LDOS, if the 
tunneling amplitude is small. When interpreting tunneling 
data measured at real tunnel junctions, one has to wonder 
which one of these two limits is the best theoretical stand- 



point. Both planar and point-contact limits are appealing 
theoretically, as they generally do not require a microscopic 
modeling of the junction. They are not equivalent, though. 
The use of planar-junction models"^'^''^"'''^ for interpreting 
STM data may therefore appear questionable, since the STM 
junction is clearly closer to the point contact than to the 
planar-junction limit. "^'^^ Planar-junction models may not 
even be appropriate for break junctions, in which the uncon- 
trolled roughness of the interface can prevent conservation 
of in-plane momentum. The similitude of the tip-sample 
junction with a point-contact is supported by first-principle 
investigations, indicating that the STM is most likely a mi- 
croscopic atom-to-atom contact.'^'^ Once formulated in terms 
of localized orbitals (e.g., Wannier functions), the tunneling 
problem may reduce in this case to a point contact involving 
only two localized orbitals, and thereby lead to a tunnel- 
ing conductance given by Eq. (6), i.e., proportional to the 
sample LDOS in the appropriate limits. 

The applications considered in the present study are sim- 
ple noninteracting single-band tight-binding models, but 
the formalism we have used has a broader validity. For in- 
stance, Eqs. (2)-(9) remain unchanged if the electrodes are 
multiband systems with interactions — the additional com- 
plexity being embodied in the real-space Green's functions — 
provided that there is no interaction between the electrodes. 
It is also possible to extend the formalism to continuum 
models; here, the main difficulty lies in the definition of 
the tunneling amplitude in Eq. (1) for continuous spatial 
coordinates. The case of superconducting electrodes re- 
quires further work. Due to the nonconservation of the 
particle number in the electrodes, additional terms involv- 
ing anomalous propagators appear in the expression of the 
current, corresponding to Andreev and Josephson contri- 
butions. This problem has been investigated for simple 
geometries using the tunneling-Hamiltonian method and 
the Keldysh nonequilibrium formalism,^'^"'^^ but the role 
of the junction' geometry for tunneling contacts involving 
superconductors remains largely unexplored. 

Nevertheless, thanks to its sub-meV resolution, tunnel- 
ing is particularly useful for investigating the excitation 
spectrum of superconductors, whose pairing gap can be 
small — typically 1 meV for a T^. of 7 K. This was beautifully 
illustrated by Giaever and his demonstration of the super- 
conducting gap,'^® and later by the detailed observation 
of vortex cores-^^"^" and impurity resonances^ ^"^"^ by STM. 
The early theory^^ has lead to the idea that the differential 
conductance is proportional to an effective DOS carrying 
information about the superconducting gap function, but 
no information about the band structure. Looking closer, 
it turns out that the information on the underlying band 
structure was discarded from the start, rather than proven 
to be absent from the tunneling spectrum. This theory also 
neglects any momentum dependence of the matrix element. 
On the contrary, the theory of Tersoff and Hamann^^ for the 
STM is based on an explicit calculation of the matrix ele- 
ment, and leads to the idea that the conductance measures 
the sample LDOS at the position of the tip. In this scenario, 
the experimental STM tunneling conductance should be 
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compared with the full LDOS containing all band-structure 
and many-body effects, as was done in recent years for 
optimally-doped cuprate superconductors of the bismuth 
family. ^"^"^^^ It was pointed out in Ref 59 that the analysis 
of STM data for quasi-2D cuprates by means of the effective 
DOS concept can lead to erroneous conclusions. 

In the Bi-based cuprates, the STM spectra develop with un- 
derdoping a more and more asymmetric shape, with excess 
weight at negative bias (occupied states).® This has been 
attributed to the onset of correlations when approaching the 
Mott insulating phase. ''""^^ Another interpretation, based on 
atomistic modeling of the tip-sample system, was recently 
proposed, in terms of copper d^i states and material-specific 
selection rules. ''•^ The asymmetric conductances in Fig. 11 
suggest yet another interpretation, based on the momentum 
dependence (rather than the energy dependence) of the 
tunneling matrix element, to explain at least part of the 
asymmetry in the cuprate STM spectra. A matrix element of 
the form (28) invariably favors the zone-center states, and 
leads for the model of Fig. 10 to an asymmetric conductance 
with excess weight at positive bias (zone-center states have 
positive energy in Sec. V because > 0). In the cuprates, 
the zone-center states lie at negative energy, so that this 
same mechanism would produce an asymmetry with ex- 
cess weight at negative bias, as observed experimentally. In 
addition to overweighting the zone-center states, the ma- 
trix element of Eq. (28) introduces local variations of this 
weighting, yielding more asymmetric conductance spectra 
in between lattice sites. If the matrix element indeed con- 
tributes to the conductance asymmetry in STM junctions to 
the cuprates, one expects to see variations of this asymmetry 
when moving the tip over subunit cell distances. Interest- 
ingly, recent investigations of the conductance asymmetry 
have unveiled such subnanometer variations, with enhanced 
asymmetry at the center of the plaquette formed by copper 
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Sites. 



VII. CONCLUSION 

We have explored the relationship between the differen- 
tial tunneling conductance and the local density of states 
(LDOS) in tunnel junctions, in particular the role of the 
junction' geometry. Our approach has been to solve exactly 
within the Keldysh formalism simple tight-binding models 
in which the geometry can be tuned from the limit of a 
point contact to the limit of a planar contact. While the 
conductance is simply proportional to the sample LDOS 
at tunneling point contacts, there can be significant differ- 
ences between the conductance and the sample LDOS at 
planar junctions, including the complete absence in the con- 
ductance of divergences present in the LDOS, which can 
be suppressed by kinematic effects (in-plane momentum 
conservation) or by interference between several tunneling 
paths. We have also studied a simple model of STM junction, 
and shown how the dispersion of electrons from the STM tip 
can induce strong variations of the local /(V) characteristics, 
and generate a contrast in the constant-current topographic 



image that is not related to variations of the LDOS. 
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Appendix A: Calculation of the current 

The tunnel current can be evaluated in the usual way,^ 
as the rate of change of the particle number in R: 



e{NR) 



(Al) 



with A/r = 2_,r ^ ' ('')''/'('')• Our convention is e = |e|, so that 
the current is positive when electrons flow from L to R, a 
positive bias V being applied to R. This corresponds to the 
usual convention adopted in tunneling spectroscopy. Since 
only J^j can change iV^, we have 

ihNi, = [N^,J^] = [NR,Ji?T] 



■S 



ni, r) [i^\l)ij(ir) - ip\r)xl,[l)] (A2) 



and the current becomes 



ir 

-2e 



(A3) 



Ir 



The quantity (\p''Xl)ip[r)) corresponds to the lesser Green's 
function at time t = 0"*".^^ It can be expressed in terms 
of the usual retarded Green's function G'^[r,l,t) = 
-i0(t)([i/)(r, t),i/)'''(I,0)] + ) and the Keldysh function 
G\r,l,t) = -i([^p[r,t),^p\l,0)]): 

{ij\l)xP[r))='-[G+[r,l,0+)-G''[r,l,0+)] . (A4) 

Introducing this into Eq. (A3), we see that, considering the 
quantities T(l,r) and G(r,[,0"'") as matrices in the {l,r} 
space, the r sum is equivalent to a matrix product and the I 
sum to a trace. Noting further that G(0+) = J G(w), we 
obtain 



e 



dwReTrr[G^(w)- G+(w)]. (A5) 



On general grounds, we have that j dcoReG'^^co) = 0. This 
can be deduced from the representation of the retarded 
function in terms of the spectral function p(e): 



ReG+(co)=a»de 



CO — e 



(A6) 
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Hence we find Eq. (2) for the current. 

In order to evaluate G^(£o) in Eq. (2), we solve the 
Dyson's equation in Keldysh space, using the same represen- 
tation as in Ref. 26: 



G+ 

Q- 



(A7) 



with G~ the usual advanced function. The Dyson equation 
for G = G(r,I,£o) expresses all ways of propagating an 
electron from r in R to f in L and reads 



G = GrTG, + GrTG, TGdTG, + . . . 



(A8) 



We have introduced the propagators of the isolated elec- 
trodes, G^, and Gj;, as well as a matrix representing the 
tunneling term, which is diagonal in the Keldysh indices. 



T 
T 



(A9) 



Equation (A8) implies matrix products in both the Keldysh 
and site indices. The term of first order in T, for instance, 
should be understood as 

^G^Cr, ri, w)T(ri, [jG^CZi, I, co), 

since the propagators G^ and G^ are restricted to the sub- 
systems R and L, respectively. Summing the infinite series 
of terms in Eq. (A8) leads to 

G = GrTG^CU - TGrTGJ"^ 

= [(GrTG^)"^ - t] ~\ (AlO) 

The matrix inversions in the Keldysh indices can be done 
explicitly, yielding the following expression for the Keldysh 
component G^: 

G^ = (n - G+TG+TT^iG+TGf + G^TG^) 

xliSl-TG-TG-T'^. (All) 

At this point, we use the assumption that each electrode 
is in thermal equilibrium and is characterized by its own 
chemical potential, say /x^ = fjLf^ + eV and pL^. The potential 
drop eV takes place in the barrier and does not affect L 
and R. Hence G^f should be evaluated for the isolated 
L and R systems at thermal equilibrium. In conditions of 
thermal equilibrium, all Green's functions can be expressed 
in terms of only two spectral functions: G^ is therefore 
not independent of G"*" and G~. Specifically, we have the 
property Gf ^(w) = [1 - 2h_i,ia))] [G+^(w) - G-^(w)] with 
R the Fermi distribution measured with respect to fx^ i^, 
and we thus obtain Eq. (3). 



Appendix B: Equivalence with Todorov et al. 

Here we show that Eqs. (2) and (3) give the same current 
as the formula obtained by Todorov et al. in Ref 22. These 



authors expressed the current in terms of the t-matrix, which 
in our notation reads 

t = T -I- TG+TG+T -I- TG+TG+rG+rG+r -I- . . . 

= r(a - G+TG+rr^ (bi) 

We first note the identity 

TrRe[-i-TG+TG+T^TG+TliG+-G'[Xl-TG-TGlT^ = 

-TrReiT\.-TG^TG+T^TiG+-G-)TG'[li-a.-TG-TGlT^. 

(B2) 

This can be proven by taking the conjugate to develop the 
real part, and then by checking order by order in T. The 
relation (B2), together with the cyclic property of the trace, 
allow to rewrite the current Eqs. (2) and (3) as 



dwTrRe(G7 - G+X^l - TG^TG^)"^ 

x(a-rG+rG+r^rG+r(2/i-2/j. (B3) 



Conjugating to evaluate the real part, and using the rela- 
tions (a - TG+TG+)-^TG+T = TG+TCl - G+rG+r)-^ 
and TG-TCn - GlTG'Ty^ = (1 - TG-TGlTy^TG'T, 
this can be rearranged as 



dwTr(G;; - G+Xa - TG~TG~r^T 
X (G+ - G-)r(II - G+TG+TrKh-fRl (B4) 



One recognizes the t-matrix of (Bl) as well as the spectral 
functions G^^ — G^^ = 2nipi j^, following the notation of 
Ref. 22, which gives Eq. (37) of Ref 22: 



/ = 



2716 



(B5) 



Appendix C: T-shaped junction for weak tunneling 

For the junction shown in Fig. 2, at lowest order in t and 
zero temperature, the conductance o"(V) of Eq. (5) reduces 
to 

„2 "i "l 



X [G-(Z2,Zi,eV)-G+(Z2,Zi,eV)]. (CI) 



Using the property L(x ± lO"*") = x =F iRe\/l — valid for 
|x| < 1, we find from Eq. (14) that 

Gl[l2,h,eV)-G+(,l2,h,eV) = 

4i/|tLl v-i 

PY2_,siii('Z'2)sin(qZi)|smq|. (C2) 
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Furthermore, we have from Eq. (11) that 

ReiG+ih - h, eV) = -Im G+(Zi - Zj, eV) ■- 
' dk 



— e'^^l^'i-'^^SCe V -2tR cos fc). (C3) 
2n 



The expression of IM^^^I^ can be worked out analytically for 
arbitrary n^; the result is displayed in Eq. (19). Finally, the 
formula (18) results by noting that 



Collecting the terms, we find 



a(V)- 



Rising I 



where we have defined 



dk 

271 



|M,j25(eV-2tRCOsfc) 
(C4) 



5(e V - 2tR cos fc) = 2TiNp{eV)5{k - k^) (C6) 



n,+l 



2sin(qZ)e''' 



(C5) with cos fco = eV/(2tjj). 
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